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Abstract 

Following the quantisation of a graph with the Dirac operator (spin-1/2) we ex- 
plain how additional weights in the spectral form factor K(t) due to spin propagation 
around orbits produce higher order terms in the small-r asymptotics in agreement 
with symplectic random matrix ensembles. We determine conditions on the group of 
spin rotations sufficient to generate CSE statistics. 
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1 Introduction 



Overwhelming evidence shows that correlations in discrete energy spectra of classically 
chaotic quantum systems generally follow the conjecture of Bohigas, Giannoni, and Schmit 
[3]. According to this the spectral statistics can be described by random matrix the- 
ory (RMT), the universality classes being completely determined by symmetry. Berry j^j 
analysed spectral two-point correlations measured with the form factor K(t) obtaining the 
leading order in its small-r asymptotics in agreement with RMT. Recent developments sug- 
gest further asymptotic terms can be attributed to correlations between classical periodic 
orbits with (almost) self- intersections [T2*l lllj. Implicit to this method is the assumption 
that the r m -term in the form factor expansion is determined only by correlations between 
pairs of orbits where the order of sections of the orbit has been changed at m — 1 (almost) 
self-intersections. The original results of Sieber and Richter, concerning systems on sur- 
faces of constant negative curvature, were extended to a class of graphs quantised with 
the Schrodinger operator by Berkolaiko, Schanz, and Whitney [21 E]- They find agreement 
with the form factor Kqq^{t) of the circular orthogonal RMT-ensemble (COE) to third 
order in r by considering correlations between orbits with up to two self-intersections. 

In jS] we showed that a graph quantised with the Dirac operator and possessing time- 
reversal symmetry produces level statistics which agree numerically with those of the Gaus- 
sian symplectic ensemble (GSE), in accordance with The form factor differs from that 
of the usual Schrodinger quantisation via weights determined by the spin transformations 
around classical periodic orbits. In the diagonal approximation we found an agreement of 
the form factor with Kgse{t) in first order. Here we show that the spin weights generate 
the additional pre-factors relating the expansions of -Kcse(t) and Kcoe(t) in all orders. 
Our approach avoids the various ad hoc assumptions introduced in the related investiga- 
tion by Heusler P and allows us to state precise conditions on the spin dynamics. The 
result also provides strong evidence for the hypothesis that the r m -term derives only from 
pairs of orbits differing in the order of sections at m — 1 self-intersections. 

2 The form factor 

In our previous work 5] we studied correlations in spectra of Dirac operators on graphs. 
Here we rather adopt the closely related point of view taken in Q] and consider the form 
factor derived from the spectrum of the S-matrix that was introduced in On the side 
of RMT we hence have to consider the circular instead of the Gaussian ensembles. 

We recall that a (compact) graph consists of V vertices connected by B bonds. The 
valency of a vertex i is v «. Let (ij) label a transition from vertex i to vertex j along a bond 
{ij}- According to [Sj for fixed energy there exist two linearly independent eigenspinors of 
the Dirac operator for each transition. The S'-matrix § is the matrix of transition elements 
connecting the eigenspinors. It is therefore a square matrix of dimension 4B. We divide 
the S'-matrix into 2x2 blocks E>w> ( kl > defining transitions between the pair of eigenspinors 
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traveling from k to / and a pair traveling from i to j. The S'-matrix is then defined by 

S«)(«> := 5 U a mki) e*<«> , (2.1) 

where is an element of SU(2) describing the spin transformation at the vertex i 

and the terms o'ajVki) define a 2vi x 2vi unitary matrix £W. The Kronecker-delta in (j2.1j) 
ensures transitions only occur between bonds connected at a vertex. The phases <f){ki} are 
random variables uniformly distributed in [0, 2ir], They define an ensemble of matrices 
§0 over which we average. Such an average is equivalent to a spectral averaging when 
considering the level statistics of a Dirac operator on the graph. 
Time-reversal invariance requires £w to be symmetric and 

§('*)(?») = ^§(*i)(fc0 j -1 _ (2.2) 

To satisfy (|2.2jl we define spin transformations 



KT • (2-3) 



The Vi elements Uj G SU(2) define all spin transformations at the vertex i. See jH| for 
details. 

Having defined the S'-matrix of a Dirac graph the form factor may be introduced as in 
[7j. We remove Kramers' degeneracy present in systems with half-integer spin and time- 
reversal invariance as in |Sj. This leaves us with iV = 2B eigenvalues of S, leading to 

The trace of § n may be expanded as a sum over the set P n of periodic orbits of length n, 

tr §™ = — A p tr K) e ^ > ( 2 - 5 ) 
peP n Tp 

where the periodic orbit p consists of a series of transitions (b\, &2> ■ ■ ■ , b n ) and 

A p e 17Tflp := cr 6nbn „ 1 cr bn _ lbn _ 2 . . . a b2bl , 

d p ■= U bnbn-l u b n -lb n - 2 _ _ u b2bl ^ 

n 
J'=l 



(2.6) 



The phases \l v are such that A p > 0, and r p is the repetition number of p so that n/r p is the 
number of possible starting positions of an orbit up to cyclic permutations. Substituting 
the periodic orbit expansion (|2.5jl into (|2.4jl and carrying out the average over we obtain 

^ympl(^) = TT^V £ tr( ^ } tr( ^ } (2 - 7) 



where r = n/2B. We label the form factor K symp i according to the symplectic symmetry 
introduced by time-reversal invariance in a system with spin- 1/2. The Kronecker-delta 
fixes contributing terms in i^ sym pi to pairs of orbits in which each bond is visited the same 
number of times. On a metric graph with rationally independent bond lengths this is 
equivalent to requiring the lengths of p and q be equal. 

For comparison the form factor studied in [21 E] for a graph quantised with the Schrodinger 
operator (spin-0) in a time-reversal symmetric fashion is 

KoMr) = ^ £ ^e*<*-*)^, (2.8) 
P,geP n p q 

where the definition of A p remains the same. This form factor is labeled by the orthogonal 
symmetry of the system. 

It was pointed out in [9j that comparing the RMT form factors of the CSE and COE 
makes clear the close connection between them, 

. . r r 2 r 3 r 4 
^cse(t) = + + + +... , 

1 2 \ 8 3 12 (2.9) 



2 -^co E y --- T + y-- + ... • 
Calling K m the term containing r m the relationship may be written 

^ \ m+l 



^cseM = ) ^coe(^) • (2-10) 

According to the conjecture of Bohigas, Giannoni, and Schmit 4j in the semiclassical limit 
we expect the form factors of quantum graphs to correspond to those of random matrices. 
In particular, 

■y \ m+l 



^symplM = ( "J ) K ^(r) . (2.11) 



It is this relation we wish to demonstrate in quantum graphs. 



3 Spin contributions to the form factor 

The form factor on graphs may be studied analytically in the semiclassical limit. The 
system is defined by the matrix § of dimension 4B and so the semiclassical limit is B — > oo. 
For small but finite r = n/2B the limit of long orbits, n — > oo, is also required. For details 
see |21 Q]. In this limit the proportion of orbits p with r p ^ 1 tends to zero so these orbits 
can effectively be ignored in equations ()2.7|) and ()2.8|) . Following pQ the sum over orbit 
pairs is organised in terms of diagrams. A diagram consists of all pairs of orbits related 
by the same pattern of permutations of arcs between self-intersections and time-reversal 
of arcs. Consequently such pairs of orbits have identical phases (fi p = <fi q . Figures [T] and El 
provide examples of diagrams. 
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The contribution to the form factor from a specific diagram D with m—1 self-intersections 

is 

KUt)--=J^ E A p A q e^-^tr(d p )tr(d q ) . (3.1) 
1 ' ( P ,q)eD n 

Here D n is the set of pairs of orbits (p, q) of length n contained in D. 

To separate spin contributions from this sum we assume that the elements d p are chosen 
randomly (independent of p) from a (sub-) group T C SU(2). This can be achieved by 
selecting the elements u? randomly from V. Then 

KyU^ = \[\^l E to K)trK)]x(^ A p A g ^~»A . (3.2) 

\J (p,q)£D n J \ (p,q)£D n J 

The second term is the equivalent contribution to the orthogonal form factor ()2.8|) . 

K£(r):=^ E AA^"^- (3.3) 

We will show that in the semiclassical limit 

1 / l\ m_1 

(tr(d p )tr(d q )r> D -.= — tr(dp) tr(^) (-- , (3.4) 

independent of the diagram D for fixed m. Substituting into ()3.2j) generates the relation 
(j2.11)) between the orthogonal and symplectic form factors. At this point we remark that 
the spin contribution to the form factor defined in (|3.4jl is different from that in [Oj. 

To determine the spin contributions we first take the case T = SU(2), ie random spin 
rotations are chosen from the whole of SU(2) with Haar measure. Each arc of p contributes 
a random element of SU(2) to d p . If (p, q) G D in d q the order of the product is changed 
and some elements are replaced with their inverse. In the semiclassical limit where the 
number of orbits tends to infinity the sum over pairs of orbits may be replaced by integrals 
over SU(2) for each arc of the diagram. To evaluate the spin contributions we require three 
identities: 

tr (xuyu) du = — - tr(xy _1 ) (3.5) 



/ 

</SU(2) 



SU(2) 



tr(xuyu x ) du = - tr(x) tr(y) (3.6) 



tr(xn) tr(yu) du = -tr(xy x ) (3.7) 

'SU(2) 2 

where u,x,y G SU(2). Equations ()3.5|) - ()3.7|) may be evaluated directly by parameterising 
SU(2). To establish ()3.4)1 we consider changing the order of arcs at one self-intersection 
or between a pair of intersections. All diagrams can be generated via these operations. 
Counting the number of intersections at which the order of elements is changed then 
proves the result. 



5 



3.1 Reordering at a single intersection 

Figure Q shows two alternative orders at a single intersection. For a pair of orbits with 
such a self-intersection the spin contribution in the semiclassical limit is 

(tr(d p ) tr(d,)> = I ... I tr(a/3Zi7<W 2 ) ^{a-f- 1 ^ 1 (3~ 1 5h) da dp d 7 d5 . . . (3.8) 

JSV(2) JSV(2) 

Changing variables so x := a(3 and y := 7<5 and using (|3.5J1 we obtain 

(tr(dp) tr(dg)) = -- / .../ ti(xkyl 2 ) tr(xl 3 yU) dx dy . . . (3.9) 

^ </SU(2) JSU(2) 

Reordering terms in d q at one vertex therefore introduces a factor of —1/2 to the spin 
contribution. 




Figure 1: A pair of orbits (p, q) with the order of arcs changed at a single self- intersection 



3.2 Reordering at a pair of intersections 

The procedure in l3.1l uses the time-reversal invariance of the system to reverse the directions 
of arcs. There is a second type of reordering of arcs independent of this symmetry which is 
possible when multiple arcs run between a pair of self- intersections, figure El The sections 
Pi . . . a 2 and 8\ . . . 72 can then be taken in either order. The relevant spin contribution is 



(tr(dp) tr(d g )) = / ... / tr(aiPiha 2 p 2 l 2 j 1 5il 3 'y 2 52h) teiai&xh^PihlxPihoiz^h) x 

isU(2) </SU(2) 

x dai da2 dPi dp 2 d.71 d72 ddi d5 2 .... 

(3.10) 

We notice that exchanging the order of the central arcs changes the order of elements of 
SU(2) at both self-intersections. From (|3.fij) and ()3.7j) it can be shown that 

/ / tr (uav ~ x bu~ x cvd) du dv = - ti(cbad) . (3-H) 

JSU(2) JSU(2) 4 

To apply this to ()3.10|) make substitutions Xj := atjPj, yj := JjSj, then 

(tv(d p ) tr(d g )) = - I ... [ tr(x 1 l 1 x 2 l 2 yihy2U) tx(x 1 l s x 2 kyihy2h) x 

4 </SU(2) </SU(2) (3.12) 

x dxi da; 2 dy± dy 2 .... 
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Reordering by exchanging arcs of a diagram introduces a factor of 1/4 to the spin contri- 
bution and requires changing the order at two self-intersections simultaneously. 




Figure 2: Reordering arcs between a pair of self- intersections 



3.3 Counting self-intersections 

All diagrams can be constructed via combinations of the procedures described in sections 
13.11 and 13.21 This is not obvious, for example a system of loops in which the same self- 
intersection is visited twice does not appear to fall within this classification, figure El In 
fact degenerate self-intersections allow both types of reordering at the intersection. To 
distinguish the cases it is necessary to follow the orbit counting each intersection when it 
is reached after determining if the order of arcs at the intersection has been changed. The 
number of self-intersections for a given diagram is then m — 1 (note our multiple counting 
of degenerate self- intersections differs from the definition in Ij). As each intersection 
effectively contributes a factor —1/2 we obtain 



(tr(d p ) tr(d,)> 




(3.13) 



The final integral over SU(2) is 

/ (tr(w)) 2 dw = l (3.14) 

isU(2) 

as the defining representation of SU(2) is naturally irreducible. 

We remark that if diagrams with m' - 1 ^ m - 1 self-intersections contributed to 
-^orth( r ) i n sucn a wa y that nevertheless K orth (r) = K coe (t) the spin contribution would 
lead to K symp i(r) ^ K C se( t )- This observation supports the hypothesis that the r m -term 
derives only from diagrams with m — 1 self-intersections. 
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Figure 3: Orbit with degenerate self-intersection 
3.4 Spin rotations from subgroups of SU(2) 

If instead of the whole of SU(2) spin transformations are chosen from a subgroup T it is still 
possible to find the connection between CSE and COE statistics. Rather than averaging 
over SU(2) the identities (J3.5J) - ()3.7|) must be understood in terms of an average over T, 
ie 




when T is finite. Let T C SU(2), viewed as a representation, be irreducible. The iden- 
tities ()3.6|) and (|3.7|) can then be derived from Schur's lemma. If and only if T is also a 
quaternionic representation (any representation equivalent to its complex conjugate but 
inequivalent to any real representation) the identity (|3.5jl may also be derived, see [Bj. As 
long as ()3.5|) - ()3.7|) hold the argument is unaffected by the use of a subgroup of spin 
transformations. CSE statistics hence depend on the subgroup of spin transformations 
providing an irreducible quaternionic representation. 

An example of a finite group of spin transformations are Hamilton's quaternions 

r = {±1, ±ia x , ±ia y , ±icr z } , (3.15) 

where <jj is a Pauli matrix. In spin transformations from this subgroup are applied to 
the cat map and CSE statistics observed. As T is both irreducible and quaternionic CSE 
statistics can indeed be expected with spin transformations taken even from such a small 
subgroup of SU(2). 

The conditions ()3.5j) - ()3.7|) depend only on the representation being irreducible and 
quaternionic consequently the argument also generalises to higher dimensional representa- 
tions of SU(2), ie higher spins. Let T C SU(2) and R S (T) an irreducible representation of 
dimension 2s + 1, ie a spin s representation, then 

^£tr(x(i?*(,)) 2 ) = 2^ytr(X) . (3.16) 

Here c = 1 for real representations and c = — 1 for quaternionic representations, however 
only quaternionic representations can have even dimension [Bj. Therefore c = — 1, compare 
(|3.5|) . implies s half-integer. ()3.6|) and ()3.7|) generalise similarly. The RMT relation between 
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the symplectic and orthogonal form factors (j2.11|) can hence be derived for half-integer spin 
provided spin-transitions generate a quaternionic irreducible representation of T. 
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